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ABSTRACT 

We present B j, R^gp a-nd /jv photometry, and optical and near-infrared spectroscopy, 
of a new ultracool white dwarf (UCWD) discovered in the SuperCOSMOS Sky Survey. 
The spectrum of SSSJ1556-0806 shows strong flux suppression due to the presence of 
collisionally induced absorption by molecular hydrogen (H2CIA), a feature character- 
istic of the cool, high density environments found in the atmospheres of ultracool white 
dwarfs. SSSJ1556-0806 therefore joins a list of <10 ultracool white dwarfs displaying 
extreme flux suppression. Synthetic model fitting suggests an effective temperature 
<3000K, which if true would make this one of the coolest white dwarfs currently 
known. We also exploit the similarity between the SEDs of SSSJ1556-0806 and the 
well-studied UCWD LHS 3250 to aid in the determination of the atmospheric param- 
eters in a regime where models consistently fail to reproduce observations. SSSJ1556- 
0806 is relatively bright {R ^ 17.8), making it particularly amenable to follow up 
observations to obtain trigonometric parallax and IR photometry. 
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1 INTRODUCTION 

The Galactic disc white dwarf population is of importance 
in obtaining disc age estimates through the use of the lumi- 
nosity function. Now devoid of the nuclear fuel that sustains 
these stars during their main sequence lifetime, white dwarfs 
shine at the expense of their residual thermal energy, slowly 
radiating away this fossil heat until they disappear from 
our telescopes as cold, compact objects. This is a very long 
process, overlapping the age of the galactic disc at the ex- 
treme. Age determination methods independent of WD cool- 
ing rates have been used to estimate the age o f the galactic 
discat around 8 Gyr (eg. I Jimenez et aL l ll998l ). ISalaris" et al.l 
l|2000l ) calculate a magnitude of M^^^j ~ 15.5 for O.GIMq 
H-rich white dwarfs that have been cooling for 8 Gyr, based 
on an interpolation of results in their Table 1. Thus, there 
is predicted to be a shortfall of white dwarfs at fainter mag- 
nitudes, a consequence of the finite age of the galactic disc 
population - white dwarfs have simply not had the time to 
cool further. The local disk white dwarf luminosity function 
is thus expected to show a marked downturn at luminosities 
corresponding to cooling ages older than that of the disk. 

Ultracool white dwarfs with effective temperatures be- 
low ~3700K fall on the cusp of this downturn in spatial den- 
sity. They represent some of the oldest stars in the sky, fos- 
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sils left over from periods of star formation in the very early 
history of the galaxy. Obtaining accurate age estimates us- 
ing this method thus depends on finding sufficient numbers 
of these extremely faint stars to constrain the observational 
determination of the luminosity function at the faint end. 
Large white dwarf surveys have in the past been hindered 
by the technical difficulties in finding such a population, 
with the limiting magnitude of the survey determining the 
cut off in the luminosity function, rather than the intrinsic 
lack of white dwarfs at cooling ages older than that of the 
galaxy. 

Modern survey class telescopes are now capable of 
observing th e faintest white dwarfs in the disc, and 
iLiebert et all l| 19881 ) were the first to detect the down- 
turn in the luminos ity function, at around MjjqJ ~ 15.76. 
IWinget et all |l983) used an early version of the observa- 
tional results to place a constraint on the age of the galactic 
disc of 9.3 ±2.0 Gyr. Currently, the largest observational de- 
termination o f the white dwarf luminosity function is that of 
[Harris et all l|2006D using SDSS photometry combined with 
USNO-B data for improved proper motions. Out of a sam- 
ple of over six thousand white dwarfs, only four objects 
fell on the downturn in the luminosity function, prevent- 
ing any quantitative conclusions as to the age of the pop- 
ulation. Evidently, the reliable identification of faint ultra- 
cool white dwarfs, even in small numbers, will enable tighter 
constraints to be placed on the spatial density of the oldest 
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white dwarfs, and the statistical uncertainty in age estimates 
to be reduced. 



Table 1. Astrometric and Photo- 
metric Data for New Ultracool 
White Dwarf SSSJ1556-0806. 



2 CURRENT OBSERVATIONAL STATUS 

In the high pressure atmospheres of cool stars colli- 
sions between H2 molecules can induce temporary dipole 
moments, allowing these sp ecies to absorb radiation 
l|Borvsow fc Frommholdlll990l ). This so-called H[2 collision- 
ally induced absorption (H2CIA) introduces broad absorp- 
tion features in the infrared, extending into the optical 
region of the spectrum in the very coolest (<3700K) at- 
mospheres. In mixed atmospheres, collisions between H2 
molecules and neutral He can induce the same process, 
and as helium atmospheres are more transparent they are 
also more dense at the same temperature, with the result 
that collisionally induced absorption becomes important at 
higher temperatures in mixed atmospheres than for pure hy- 
drogen atmospheres. The photospheric opacity of UCWDs is 
thus dominated by this process, resulting in spectra showing 
a strong flux deficit towards the red relative to a blackbody 
form. 

Cool white dwarfs can be roughly divided into two 
classes. Below around 6000K, opacity due to H2CIA is rel- 
atively mild and leads to flux suppression in the infrared 
region of the spectrum. Objects in this category include the 
white dwarfs F351-50, LHS 1126 a nd the well characte rised 
halo white dwar f WD 0346-H246 (|Hamblv et al.lll997l ). for 
which iBergeronI l|200ll ') calculates a T^q of 3780K. With 
decreasing temperature, H2CIA becomes a more dominant 
source of photospheric opacity. 

Below around 3700K, flux suppression extends into 
the optical region of the SED, shifting the peak blue- 
wards with decreasing Tgg-. This results in a reversal of 
the slow redwards trend in the optical and infrared colours 
of cooling white dwarfs as this process begins to dominate 
the atmospheric opacity. Such ultracool objects are often 
misident ified by a utomated search routines, a fact used by 
[Gates et al.l (|2004l ) to locate five white dwarfs in the Sloan 
Digital Sky Survey that show varying degrees of flux sup- 
pression. Several other UCWDs have been discovere d in the 
SDSS, including SDSS J 1337 bv iHarris et al] (|200ll ). These 
are generally very faint {r > 19), making them difficult ob- 
jects to study in terms of parallax follow up or IR photom- 
etry. 

Extreme examples of flux suppression include the stars 
LHS 1402 and LHS 3250. Obs e rvatio ns of the latter are pre- 
sented bv lOppenheimer et al.l (|200lh alongside those of WD 
0346+246, exhibiting the difference between these two rough 
groupings. 



Property 



Value 



Designation 
Right Ascension^ 
Declination^ 
Epoch 

Mtot (as yr-^) 

Position Angle 

Bj 

R59F 

In 



SSSJ1556-0806 
15 56 47.3181 
-08 5 59.713 
1992.480 
0.421±0.0055 
117.913 
18.800 
17.840 
18.058 



^ Coordinates given in equinox 2000 



covers the entire sky, and with ~5 times more sky coverage 
than SDSS provides a valuable new resource for identifying 
UCWDs. 

We have obtained high proper motion candidates (O'.'lS 
yr~^ ^ /i|^Ql^ ^ lO'.' OO yr~^) by use of an object pairing algo- 
rithm described in lHamblv et al.l (|2004l ). White dwarf can- 
didates are selected based on their Reduced Proper Motion 
(RPM), an estimate of absolute magnitude that uses proper 
motion as a distance proxy. A plot of RPM against colour 
is topologically equivalent to the classical HR diagram, and 
white dwarfs can be identified in the usual manner by their 
extreme subluminosity at a given colour. The use of the 
SSS in defining a clean, complete sample of white dwarfs 
based on the t echniq ue of RPM is described thoroughly in 
iHamblv et all l|2005l ). 



4 A NEW ULTRACOOL WHITE DWARF 

We have identified a new ultracool white dwarf in data taken 
from the SuperCOSMOS Sky Survey. This star, designated 
SSSJ1556-0806, shows the strong H2CIA fiux deficit asso- 
ciated with the cool, high pressure photospheres of white 
dwarfs of effective temperature <3700K. 



4.1 Observations 

SSSJ1556-0806 was identified as a white dwarf based on it's 
RPM, as described in section|3] It was selected as an UCWD 
candidate from it's location in the two-colour plane (see Fig- 
ure [T| , below the main locus of white dwarfs and close to 
a group of known UCWDs. Various astrometric and photo- 
metric data are presented in Table 1. 



3 THE SUPERCOSMOS SKY SURVEY 

The SuperCOSMOS Sky Survey (SSS) has been compiled 
by scanning photographic Schmidt plate images taken over 
a baseline of up to fifty years, in three different photographic 
pass bands and one twice , giving four epochs of observation 
fsee lHamblv et aLll2001al .b.c). The photographic passbands 
Bj, R^gF and 7jv differ from standard Johnson BRI, colour 
transformations are dealt with in section \5l2\ The SSS now 



4.2 Spectroscopy 

Follow up spectroscopic observations are necessary to con- 
firm the identity of SSSJ1556-h0806 as an UCWD. These 
were carried out over two semesters at two different insti- 
tutions. Optical spectroscopic observations were made using 
the Marcario Low Resolution Spectrograph on the Hobby- 
Eberly Telescope at McDonald Observatory, Texas, on the 
night of 2006 July IS^'^. Grism 2 in conjunction with a I'.'S 
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Figure 1. SSSJ1556-0806 lies below the main white dwarf locus 
in the colour-colour plane, in roughly the same region as sev- 
eral known ultracool white dwarfs that display extreme flux sup- 
pression. Open triangles represent previously confirmed UCWDs 
that appear in SSS. Open squares are normal white dwarfs that 
appear in SSS and h ave been spectroscopically confirmed by 
iMcCook fc SionI l|l999h . The filled triangle is SSSJ1556-0806. Er- 
ror bars have been omitted from the normal white dwarfs for 
clarity; they are of order aB — R,R—i ~ 0.1. 
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Figure 2. Spectra of SSSJ1556-0806 taken using both the Hobby- 
Eberly telescope and William Herschel telescope. The break in 
the WHT spectra is due to the dichroic beam-splitting mirror 
used in the dual-arm ISIS spectrograph. The flux scale is accu- 
rate for the WHT spectrum; the HET spectrum has been offset 
vertically for clarity The SDSS spectrum of UCWD LHS 3250 has 
been included to demonstrate the similarity between the SEDs of 
these two objects. 



vious models. Only pure hydrogen (qjj 



10" 



and pure 



helium (qjj = 10 ) atmospheres are considered, mixed 
atmosphere models have not been used here. 



slit produced spectra with a resolution of ~6A over the range 
4280- 7340 A. Optical/near-IR observations were made at the 
William Herschel Telescope in the Canary Islands on 2007 
April 24*^^ using the dual-arm ISIS spectrograph. A l'.'5 slit 
was used together with the R300 grating on the blue arm and 
R158 on the red arm. The resolving power of each arm was, 
respectively, ~9A and ~llA over the ranges 3160-5270A 
and 5710-10350A. Data reductions were carried out using 
standard routines in IRAF. In Figure [2] we present the spec- 
tra obtained at both of these observing runs, offset vertically 
for clarity. The SDSS spectrum of the well-studied UCWD 
LHS 3250 has been included (offset vertically) to show the 
similarity between the spectral energy distributions of these 
objects, a feature the will be used to aid analysis in section 
[531 

5 ATMOSPHERIC PARAMETERS 
5.1 The models 

In order to derive rough atmospheric parameters for this 
star, we have fitted the photometric colours to grids of 
DA and DB white dwarf sy nthetic colours computed by P. 
Berge ron and described in [Fontaine. Brassard fc BergeronI 
l|200ll ). The DA grid ranges in effective temperature from 
1500K to 100,000K, and the DB grid from 3500K to 
30,000K; both cover a range of \og{g) from 7.0 to 9.0. A va- 
riety of models at fixed masses have also been fitted, these 
range in mass from 0.2Mq to 1.2Mq in steps of O.IM© 
and cover a range of Tgg- and log{g) similar to the pre- 



5.2 Colour transformations 

The magnitude system used in the SSS is non-standard. The 
subscripts attached to the magnitudes Bj, R^gp and In 
refer to particular combinations of colour filter and photo- 
graphic emulsion, resulting in response functions different 
to the standard Johnson BRI. In particular, Bj deviates 
significantly from Johnson B, being closer to a Johnson V , 
and i?59F lacks the long red tail that is a characteristic of 
Johnson R. 

In order to fit the photographic magnitudes of 
SSSJ1556-0806 to the photoelectric magnitudes of the mod- 
els, transformation of UKST photographic magnitudes was 
therefor e necessary. For this we used the colou r equati o ns de - 
rived bv lSalim et al] (|2004l ') for B and /, and lBesselll (Il986l ) 
for R. Note tha t othe r colour transformations exist, e.g. 
iBlair fc Gilinor^ (Il982l ) derive different transformations us- 
ing all standard stars in a given region, however Salim et al. 
use a sample of cool white dwarfs and hence obtain equa- 
tions more suited to our purpose. 

5.3 Fitting procedure 

We have transformed the photographic BjR^gplN magni- 
tudes of SSSJ1556-0806 to Johnson magnitudes as described 
previously, and used a variance weighted least squares 
method to fit these to the photometric colours of the models. 
The model predictions are provided in terms of broadband 
fluxes, hence the fits are done using these alone rather than 
the full spectrophotometry. 



4 N. R. Rowell, M. Kilic and N. C. Hambly 




4000 5000 6000 7000 8000 9000 
Wavelength (A) 

Figure 3. Photometric fluxes of the best fltting DA (circles) and 
DB (diamonds) models are plotted over those of SSSJ1556-0806 
(boxes). Large error bars are due to a combination of absolute 
photometric uncertainty and uncertainty in transformations to 
Johnson passbands. The dashed line is a ninth order polynomial 
fit to the continuum flux of SSSJ1556-0806. Models are provided 
in terms of broadband fluxes, so our fitting procedure is on these 
alone rather than the full spectrophotometry. Model fluxes have 
been normalised to the R band at 6OOOA. The solid line is a 
blackbody fit to the blue end (3500-4800A) of the SSSJ1556- 
0806 spectrum, the region of the optical spectrum of UCWDs 
that most closely approximates a blackbody, and is intended to 
be roughly illustrative of the effective temperature range. 



As the atmosphere content is uncertain, we have found 
the best fitting DA and DB models for comparison. The 
mass distribution of white dwarfs is tightly distributed about 
the mean, and so for the DA and DB fits we have restricted 
the mass to M = 0.603±0.081 Mo and M = O . nSzfc O.lllMa 
respectively, values taken from lKepler et all (|2007l ). 

The parameters of the best fitting models are given in 
Table 2, and a fit of the model colours to the stellar colours is 
presented in Figure [S] with magnitudes normalised to the R 
band. These have been converted to fluxes, allowing visual 
comparison with the WHT spectrum, using the prescrip- 
tion outlined below. Magnitudes in different passbands have 
been converted into approximate average fluxes at central 
wavelengths using the equation: 

m = -2.51og/r + C^ (1) 
where 



5.4 Temperature and atmosphere content 

For the DA model, the fltted effective temperature of 2500K 
agrees qualitatively with what one would expect based on 
visual examination of the SED - namely, that at this tem- 
perature one would expect strong flux suppression extending 
into the optical region of the spectrum, and a lack of any 
discernible atomic absortion lines. The temperature of the 
DB model is suspiciously high - even though H2CIA sets 
in at higher temperatures in mixed atmospheres, at 7000K 
one would not expect to see such strong fiux suppression. A 
cooler temperature is also suggested by fitting a blackbody 
curve to the blue end - the region of the optical spectrum in 
UCWDs that most closely approximates a blackbody form. 
Our least-squares fit (in steps of 50K in the range 3500- 
4800A) implies an effective temperature <4000K, at odds 
to that predicted by the DB model. Note that the fitting 
of a blackbody is not in general an accurate way to mea- 
sure stellar temperature, however in a temperature regime 
where stellar models also fail to accurately reproduce obser- 
vations (see section 15. 5p this can at least provide an illus- 
trative temperature measurement, and is used frequently in 
the literature. 

The age of the DA model would make this star several 
Gyr older than the galactic disc, and thus a strong candidate 
for a member of the halo population. However, the calcu- 
lated tangential velocity is not consistent with halo kinemat- 
ics, and would suggest membership of the disc population 
instead. The DB model places the star at sufficient distance 
that the tangential velocity is indicative of halo membership, 
and is comparable to that of the classic halo white dwarf 
WD 0346-1-246 which has vt ~ 170 kms'^ While the cool- 
ing age of ~2.3 Gyr is considerably younger than the halo 
population, we cannot rule this out as a halo star based on 
this alone. The total stellar age (MS -f WD) is required, 
and this can only be determined with accurate distance and 
mass measurements. 

One expects UCWDs to have hydrog en domin a ted at - 
mospheres due to arguments outlined in iBereeronl l|200lh . 
namely that during their long cooling periods, these stars 
travel throughout the galaxy accreting material from the 
interstellar medium. While heavier elements tend to quickly 
sink below the photosphere, accreted hydrogen will float on 
top, resulting in an atmospheric hydrogen abundance that 
can only increase with time. Combined with the observa- 
tional fact that, at least for hot white dwarfs, DAs outnum- 
ber DBs by a factor of ~10, we regard the DA model as 
a more realistic solution. Indeed, the DB model provides a 
poor fit to the photometric data and can be ruled out as 
a likely solution. However, with only two colour indices to 
fit, it is impossible to draw any firm conclusions as to the 
nature of this object. 



is the average flux in band m, /(A) is the flux received from 
the star at wavelength A, and 5m (A) is the transmission 
funct ion for bandpass m. Th e constants Cm have been taken 
from [Bergeron et al.l l|l997h , and are calculated from Vega 
fluxes. 



5.5 Comparison with LHS 3250 

Synthetic WD spectra fail to accurately reproduce the SEDs 
of UCWDs, suggesting inadequate or incomplete physics 
used in the model atmosphere calculations. In particu- 
lar, pure hydrogen models show a large absorption feature 
around 7500A which is not seen in the spectrum of any con- 
flrmed UCWD, and the peak of the flux distribution is al- 
ways predicted too narrow. Note that here we have used 



Table 2. Atmospheric parameters derived for SSSJ1556- 

0806 BY FITTING DA AND DB MODEL ATMOSPHERES. 



Spectral Type 

Property DA DB 
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Table 3. Atmospheric parameters for 
SSSJ1556-0806 derived using LHS 3250 AS 

AN ANALOGUE. 

Property Value 

Effective Temperature (K) <4000K 

Atmosphere content He rich 

Ms, 16.302 ± 0.367 

Mn^gj, 15.274 ± 0.369 

Mj„ 15.730 ± 0.392 

Distance (pc) 31.6 ± 1.4 

Tangential Velocity (kms^^) 63.1 ± 2.8 



Effective Temperature (K) 


2500 


7000 


log(9) 


8.04 


8.199 


Mass 


0.6 Mq 


0.7 Mq 


Bolometric Correction 


0.364 


-0.124 


Mbol 


17.942 


13.702 


Age (Gyr) 


12.47 


2.315 


Ms 


18.905 


14.198 


Mn 


17.023 


13.572 


Mi 


17.149 


13.320 


Distance (pc) 


12.85 ± 0.15 


80.71 ± 5.52 


Tangential Velocity (kms^-"^) 


25.64 ±0.16 


161.06 ± 8.08 



only model colours, however these are calculated by inte- 
grating model spectra over the appropriate bandpasses, and 
are assumed to suffer the same problems. 

The similarity between both the SEDs (see Fig[2| and 
photometric colours {Bj - Rf = 0.96 and 1.028, Rf - In 
= -0.218 and -0.456, respectively) of SSSJ1556-0806 and 
the UCWD LHS 3250 motivates us to consider an alterna- 
tive approach to the analysis of S SSJ1556-0806 . LHS 3250 
has been studied ri g orous ly - see [Harris et al.1 ||1999D and 
[Bergeron fc LeggettI l|2002h . Bergeron and Legget fit a vari- 
ety of atmosphere models, covering a large parameter space, 
to their data, which include trigonometric parallax, JHK 
photometry and low resolution optical spectra. Their fig. 3 
demonstrates the disagreement between observed and theo- 
retical SEDs. On this basis they rule out a pure H solution 
for LHS 3250, instead finding a mixed atmosphere solution 
with log N(H) /N(He) = -4.7, log (g) = 8.00 and Tgfj =3480K. 

To proceed with our analysis of SSSJ1556-0806, we shall 
assume that these stars share similar absolute magnitudes 
and atmospheric properties. The distance to LHS 3250 is 
well constrained from trigonometric parallax measurement, 
allowing us to convert it's SSS apparent magnitudes to abso- 
lute magnitudes in the photographic BjR^gFlN passbands. 
The apparent magnitudes of SSSJ1556-0806 thus place it 
at a distance of 31.6 ± 1.4pc, on combining the distances 
calculated from the individual passbands. This suggests a 
tangential velocity of 63.1 ± 2.8kms~^, consistent with the 
galactic thick disc population. 

The atmospheric and evolutionary properties are more 
tricky; Bergeron & Legget point out that the inconsistency 
between observed and theoretical SEDs rules out precise de- 
termination of these parameters. In the circumstances, we 
can do no more than say that this star most likely has an 
effective temperature below 4000K and a helium-rich atmo- 
sphere. These various properties are summarised in Table 3. 



6 CONCLUSION 

We have presented observations of a new ultracool white 
dwarf found in the SuperCOSMOS Sky Survey, and made 
an initial effort to understand the unusual SED of this ob- 
ject using both atmosphere models for cool white dwarfs and 



comparison to the well-studied UCWD LHS 3250. Without 
better photometry, models are poorly constrained and can 
do no more than give a rough interpretation of the object. 
The current inconsistencies between the theoretical and ob- 
served spectral energy distributions of UCWDs rules out 
any precise determination of the atmospheric parameters. 
The most likely interpretation of this object is as a thick 
disk white dwarf, with a helium dominated atmosphere and 
an effective temperature <4000K. 

Infrared photometry is essential for a good fit to theo- 
retical models. Trigonometric parallax observations can also 
provide constraints by measuring both the absolute magni- 
tude and, indirectly, mass. We currently have time allocated 
on the robotic Liverpool Telescope for a program to obtain 
parallax measurements of UCWDs. SSSJ1556-0806 has been 
added to our queue-scheduled list and has already had one 
epoch of observation carried out at high parallax factor. 



REFERENCES 

Bergeron P. 2001, ApJ, 558, 369 
Bergeron P., Leggett S. K. 2002, ApJ, 580, 1070 
Bergeron P., Ruiz M. T., Leggett S. K. 1997, ApJS, 108, 
339 

BesseU M. S. 1986, PASP, 98, 1303 
Blair M., Gilmore G. 1982, PASP, 94, 742 
Borysow A., Frommhold L. 1990, ApJL, 348, L41 
Fontaine G., Brassard P., Bergeron P. 2001, PASP, 113, 
409 

Gates E., et al. 2004, ApJL, 612, L129 

Hambly N. C, et al. 2001a, MNRAS, 326, 1279 

Hambly N. C, Irwin M. J., MacGillivray H. T. 2001b, MN- 
RAS, 326, 1295 

Hambly N. C, DavenhaU A. C, Irwin M. J., MacGillivray 
H. T. 2001c, MNRAS, 326, 1315 

Hambly N. C, Digby A. P., Oppenheimer B. R. 2005, 14th 
European Workshop on White Dwarfs, 334, 113 

Hambly N. C, Henry T. J., Subasavage J. P., Brown M. A., 
Jao W.-C. 2004, AJ, 128, 437 

Hambly N. C, Smartt S. J., Hodgkin S. T. 1997, ApJL, 
489, L157 

Harris H. C, Dahn C. C, Vrba F. J., Henden A. A., Liebert 
J., Schmidt G. D., Reid I. N. 1999, ApJ, 524, 1000 

Harris H. C, et al. 2001, ApJL, 549, L109 

Harris H. C, et al. 2006, AJ, 131, 571 

Jimenez R., Flynn C, Kotoneva E. 1998, MNRAS, 299, 
515 



6 N. R. Rowell, M. Kilic and N. C. Hamhly 



Kepler S. O., Klcinman S. J., Nitta A., Koester D., Cas- 

tanheira B. G., Giovannini O., Costa A. F. M., Althaus 

L. 2007, MNRAS, 375, 1315 
Liebert J., Dahn C. C, Monet D. G. 1988, ApJ, 332, 891 
McCook G. P., Sion E. M. 1999, ApJS, 121, 1 
Oppenheimer B. R., et al. 2001, ApJ, 550, 448 
Salaris M., Garci'a-Berro E., Hernanz M., Isern J., Saumon 

D. 2000, ApJ, 544, 1036 
Salim S., Rich R. M., Hansen B. M., Koopmans L. V. E., 

Oppenheimer B. R., Blandford R. D. 2004, ApJ, 601, 1075 
Winget D. E., Hansen C. J., Liebert J., van Horn H. M., 

Fontaine G., Nather R. E., Kepler S. O., Lamb D. Q. 1987, 

ApJL, 315, L77 

This paper has been typeset from a T^iX/ I^T^iX file prepared 
by the author. 



